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1 Introduction 

In this paper we consider the following example of unbounded lattice spin sys- 
tem. Let A C 1i d be a finite set of lattice sites which contains the origin 0. For 
each site x£ Awe associate a complex- valued random variable ip{x) called a 
spin or a field. Note that we will use "*" to denote complex conjugation instead 
of a bar. The collection ip\ = (^(x)) x£ a £ C A of these variables is sampled 
according to the finite volume Gibbs measure 

(.) A = -L [ Df^ e -^(«(.) (i) 

where 

Z A = I Di/j*Dip e ~ H ^ A) > (2) 
xeA 

is proportional to the Lebesgue measure in C A , and the Hamiltonian with free 
boundary conditions is given by 

H K m= £ J(x-y)r( X )^(y) + ^J>(x)| 4 . (3) 

(x,y)eA 2 xeA 

The assumptions on the parameters appearing in the Hamiltonian are the fol- 
lowing. 

• A > 0. 

• The pair potential J is a function Z d —>■ C of compact support such that 

J(x)* = J(-x) . 
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• = Yl l^( x )l i s finite and moreover satisfies < < J(0). 

Let ip^ stand for either ip or ip*. The main goaf of this article is to study 
the so called truncated correlations functions, defined by 

(^(x 1 ),...,^(x n ))J = 
' : Log((l + a 1 ^(xi))...(l + a ri V tt (xn)))A 



da i 9a n 



ai,...,a„=0 



where Log denotes the principal logarithm of a complex number, i.e. Log(re ie ) = 
ln(r) + iO, with r > and £ (— tt, it]. It is a standard task to show that the 
series in ai,...,a n is analytic in a small poly disc around cti,...,a n = 0. 
These truncated correlations are also known in the literature as cumulants, 
semi- invariants or connected correlation functions. 

We will find uniformly in A and the assignments for the (j symbols, explicit 
bounds of the form 

£ K^ 1 (0),^ 2 (x 2 ),...,^"(x„))H <c„(J,A) (4) 
x2,...,x n eA 

also called ^-clustering estimates. We will derive such bounds in four different 
regimes: 

• J(0) big, i.e., the large-mass regime; 

• small, i.e., the regime of small interaction between lattice sites; 

• A big, or the 'high temperature' regime; 

• A small, i.e., the near-Gaussian regime, or the 'low temperature' regime. 

Remark 1 This way of introducing the distinction between high and low tem- 
perature is not what is usually done in textbook treatments of the Ginzburg- 
Landau theory of phase transitions (see e.g. [34, Ch. 2]) which addresses models 
similar to the one above. Instead, we introduce the temperature dependence by 
replacing the Hamiltonian H\(ip\) by (3H\(ip\) where (3 as usual denotes the 
inverse of the temperature times Boltzmann's constant. However, one can do 
the change of variable ip = f3~^ip' and absorb (3 into a modification of the cou- 
pling A — ► /3 _1 A. In this setting, A big coincides with high temperature, and A 
small with low temperature. In the present situation of a massive lattice model, 
one does not expect a phase transition at 'low temperature '. 

In the statement of the theorems below we will denote by the same symbol 
0(1) the various constants which appear. These are absolute constants such as 
y/n or quantities which only depend on the dimension d which is fixed through- 
out. In the following sections we will express such constants in details, to prove 
that they are effectively computable. However we are not interested in this 
paper in optimal bounds so the values that we give may be improved with more 
refined estimates. 

Our results for the various regimes above consist of the following theorems. 
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Theorem 1 In the event that > and A > are fixed, there exists a 
K(J^,X) > > which is independent of the volume A such that the l 1 - 
clustering estimate 



|(V' ttl (0),^ 2 (x 2 ),...,V tt "(x ri ))I| <n! 

x 2 ,...,x„eA 
holds as soon as J(O) > K(J^,X). 



0(1) 



VW) - J* 



Theorem 2 In the event that J(O) > and A > are fixed, there exists an 
e(J(0),A) G (0, J(O)) which is independent of the volume A snc/i i/iai i/ie Z 1 - 
clustering estimate 



(0),^ 2 (x 2 ),...,VHx n )>H <n! 



x2,...,x n eA 



/io/ds as soon as < < e(J(0), A). 



Theorem 3 in i/ie et>en£ i/iai J-l > and J(0) > are /ixed, i/iere exists 
a K(J(0), J^) > which is independent of the volume A such that the l 1 - 
clustering estimate 

\W il (0),^(x 2 ),...,^(x n ))l\<n\y. [0(l)xA-i]" 

x 2 ,...,x„eA 
holds as soon as A > K(J(0), J^). 

Theorem 4 Let N > 1 be an integer and suppose that J is fixed. Then there 
exist quantities e(J) > ; ci(N, J) > and c 2 (J) > such that, for any A with 
< A < e(J), for any nonempty finite volume A C 1? , and any even integer n, 
n > 2(N + 1), we have the ^-clustering bound 

Y |<^ 1 (0),^ 2 (x2),...,V' ttn (xn))Il <ci(N,J) x c 2 (J) n x X N x n\ . 
x2,...,x„eA 

Theorem 5 There exists a quantity cs(J) > such that for the same e(J) as 
in the previous theorem, for any A with < A < e(J), and for any nonempty 
finite volume A C 7j d , one has 

£ |(^ 1 (0),V fe (x 2 ))lA-(^ 1 (0),^ 2 (x 2 ))I,ol <c 3 (J)A 
x 2 eA 

where we have restored the X-dependence in the notation for truncated correla- 
tion functions. 

Such ^-clustering estimates in statistical mechanics have a long history. 
See [45, 46] for the case of classical gases at low activity. The importance of such 
estimates in relation to analyticity of thermodynamic functions was stressed 



4 



in [19]. However, to the best of our knowledge, clustering estimates which 
include the results of the previous theorems have not previously appeared in the 
literature. The cluster expansion methods we used to derive these results also 
provide the exponential tree decay of the truncated correlations as a function 
of the locations x l5 . . . ,x n of the sources. The reader can easily extract such 
decay estimates from the proofs provided here. However, in this article we focus 
instead on I 1 estimates where the Xj are summed over the lattice. Indeed, our 
primary motivation comes from the recent work of Lukkarinen and Spohn [36, 
37] where they obtained the rigorous control of the kinetic regime for the time- 
evolution of a class of interacting quantum fluids. Their result is conditional on 
establishing ^-clustering estimates such as the ones provided in Theorem 4 and 
Theorem 5. Similar estimates may also be obtained in the Fermionic case in 
an easier way (see [48]). Our methods are robust enough to handle much more 
general undounded spin systems, but for better readability we refrained from 
stating our results with maximal generality, and restricted our attention to the 
complex Bosonic model presented above which is the one needed for [36, 37]. 

Another motivation for the present work comes from the recent interest in 
the decay of correlations for unbounded spin systems, especially in relation to 
Log-Sobolev inequalities (see e.g. [59, 57, 6, 43, 29, 28] and references therein). 
For the equivalence between the exponential decay of the truncated 2-point 
function, the Log-Sobolev inequality and the spectral gap property, in the case 
of unbounded spins, see [58]. The study of the decay properties for the 2-point 
function for this class of models has a long history [33, 26, 31, 17, 52]. In the 
low temperature case, i.e., in a regime similar to the setting of Theorems 4 and 
5, recent proofs of exponential decay for the truncated 2-point function were 
given in [50, 7, 51, 8, 40]. These works use the methods introduced by Helffer 
and Sjostrand in [30]. However, to the best of our knowledge, one has not 
been able to treat higher truncated n-point functions with Helffer-Sjostrand 
and Witten laplacian techniques. The only such results [30, 35] for higher 
correlations that we are aware of concern centered moments, i.e., expectations 
of the form ((Xi — (Xi)) . . . (X n — (X n ))) which are different from fully truncated 
correlations (Xi, . . . , A n ) T . 

It is quite well known that the first three regimes mentioned above are 
amenable to cluster expansion techniques of the kind that is standard in the 
statistical mechanics literature (see e.g. [46, 39, 49]). Much less known, in the 
mathematical analysis and probability theory communities, is that the low tem- 
perature regime can also be treated using a special kind of cluster expansion 
technique. We refer to the latter as the field theoretic cluster expansion. It 
also has a long history, and originates in the work of Glimm, Jaffe and Spencer 
in constructive quantum field theory [23, 24]. It has then been simplified and 
improved by many authors (see in particular [12, 38, 41, 22, 44] and references 
therein). Our approach for the proof of Theorems 4 and 5 owes much to the 
pioneering work [20] in the context of P((f>)2 quantum field theories. This was 
adapted to the lattice setting in [18, 55]. These use the original Glimm- Jaffe- 
Spencer cluster expansion, which involves a decoupling procedure for the Gaus- 
sian measure. A simpler way to do this was introduced by Brydges, Battle and 
Federbush [14, 10] and it is based on a combinatorial tree expansion identity. 
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The third generation of such decoupling procedures, which is the one used in 
this article, is based on the Brydges-Kennedy-Abdesselam-Rivasseau (BKAR) 
forest formula [15, 2]. We give an account of this basic tool and of the general 
results of cluster expansion in §2.4 and §2.5. We will use it many times, in 
the proof of the estimates, first for the 'high-temperature' regimes and then for 
the 'low temperature' regime. One of the difficulties which we have to address 
here and which does not seem to have received a great deal of attention in the 
previous literature, is the aim for bounds which are uniform in n, growing as n\, 
and where one simultaneously extracts as many powers of A as possible. Such 
an extraction of perturbation theory in the context of constructive field theory 
can be done using additional Taylor expansions (see e.g. [9, §5.14]). However, 
we have been unable to extract the optimal bound C n \ 2 ~ l n\ predicted by tree 
level perturbation theory, in a way which is uniform in n. The precise statement 
of this problem is given as Conjecture 1 from §2.2. Of course, one can ask the 
same question for the real-valued scalar field model with ip A interaction, on the 
lattice. Note that a related bound with a factor A 4 was obtained by Brydges, 
Dimock and Hurd [13, Theorem 9] in the context of the UV limit of the ^3 
model. 

Let us conclude this introduction by indicating some of the notation used 
throughout this article. We use | • | for the cardinality of finite sets. If n is a 
nonnegative integer, the set {1,2,..., n} is denoted by [n]. We will denote by 
!{• • •} the characteristic function of the condition between braces. 



2 Preliminaries 

2.1 Basic properties of the model 

By hypothesis on the function J, the matrix J = ( J(x — y)) x ,yeA is Hermitian 
positive definite. This can be proved as follows. Let us denote the inner product 
on / 2 (A) by < ipi,ip2 >= SxeA V'i( x )V ; 2( x )) and the norm of a vector by 
|| • ||. The hypothesis J(x)* = J(— x) trivially implies that the matrix J = 
(J(x — y)) x ,yeA is Hermitian jt = J. Besides for any field i/j on / 2 (A), on has 

<i)Ji)> = E V>*(x)J(x-y)^(y) 
x.yeA 

= J(0)\M\ 2 + E V>*(x)J(x-y)V(y) 



x,ygA 
X5^y 



and 



E Y>*(x)J(x-yMy) 



x.ygA 
x#y 



< E (Wx)|vV(x-y)|) (Wy)lvV(x-y)|) 



x.ygA 
x^y 



<\ E IV>(x)lV(x-y)| + J E IV>(y)lV(x-y) 



2 ^ ITV " v Jn 2 

x.ygA x.ygA 
x#y x^y 



6 



where we used the inequality ab < \{a? + b 2 ). Therefore the last expression is 
bounded by J^||^|| 2 and 

<^,J^>>(J(o)-J^)M\ 2 

so that J is positive definite by the hypothesis J(0) > J^. The covariance 
matrix C = J -1 = (C(x, y)) x ,yeA is well defined, and there exists a unique 
mean zero normalized Gaussian probability measure denoted by dfici^P* >Y0 on 
C A with covariance matrix C, i.e., such that 

/ ^ c (Y^V)^(x)V*(y) = C(x,y) , 
Jc A 

I dnc{r^)r{*W{y) = o, 

[ dficir,^) ^(x)V(y) = o. 

The measure can be written 

dnc(il>*,il>) = (detJ)e-^y^ 2j{x ~ y)r{x) ^ y) DrDij . 

The moments of this measure can be expressed via the Isserlis-Wick Theo- 
rem [32, 56] as follows: 

• If p ^ q then 

/ YO ^(xi) . • • V(xp)V*(yi) • ■ ■ Y>*(y ? ) = o . 

Jc A 

• In the p = q case, one has 

/ dficir^) V>(x0 . . . ^(x p )V*(yi) . . . V*(y P ) 

= C(x 7(1) ,yi)...C(x 7(p) ,y p ) , (5) 

i.e., one sums over all possible pairwise Wick contractions of the ^'s with 
the ifr* 's, as indicated by the permutation 7. 

2.2 Feynman diagrams and tree level analysis 

The relation (5) may be expressed in terms of oriented graphs (Feynman dia- 
grams) in which for each ^>(x) we draw an ingoing half edge from the vertex x, 
for each ip* (y) we draw an outgoing half edge from the vertex y , and the expec- 
tation is the sum over all the possible oriented graphs obtained contracting only 
outgoing half edges with ingoing half edges and associating a free propagator 
C(x, y) to each edge. 

In this section we will focus on the near-Gaussian regime, i.e., A small, and 
we will analyze the truncated n-point function (ip^ 1 (xi ),..., Y>" n ( x n))J L \- A 
trivial change of variable ip — > e td il> shows that truncated correlation functions 
vanish unless the number of arguments n is even, and there are an equal numbers 
of ^'s and ip* 's involved. We will therefore always assume this to be the case. 

An easy but crucial lemma we will later need is the following. 
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Lemma 1 For n even, n > \, and for all k < ^ — 1 we have 

. 



^)VHxi),...,V^(x„)>I A 



A=0 



We assume the reader is familiar with the rigorous formalism of Feynman 
diagrams used to express formal perturbation theory. A precise mathematical 
treatment can be found in [1, 47]. A well-known fact from this formalism is that 
the function (^"(xi), . . . , ^(x„))j[ A is C°° in A in the interval [0, +oo), and that 
its Taylor series at the origin, seen as a formal power series in A, is the sum of the 
contributions of all connected Feynman diagrams with external legs x 1; . . . , x„. 
It is a trivial exercise in organic chemistry to see that the minimal number N 
of internal 4-valent vertices needed to build such a connected graph is § — 1- 
This corresponds to tree graphs. The lemma is an immediate consequence of 
this fact. 

The remainder of this section will be devoted to some heuristic consider- 
ations which we hope will shed some light on the near-Gaussian regime. It 
is part of constructive field theory folklore that one should expect the correct 
bound in the ^-clustering estimate to be dictated by the contribution of tree 
graphs. For a given tree graph, the I 1 sum over the sites X2, . . . ,x n in Z d can 
be easily bounded, using a standard pin and sum argument, by Cst n . This is 
because the hypotheses on J imply the exponential decay of the free propaga- 
tor C(x, y), as is recalled in §2.3. The issue is the number of trees. Because 
of the constraints due to the orientations of the edges, the counting is not an 
immediate consequence of Cayley's formula. One can nevertheless obtain an 
exact formula. 

For n > 2 and even, let denote the number of Wick contractions 7, 
with n external legs and N = ^ — 1 internal vertices, which produce connecting 
trees. One can easily check by inspection that k\ = 1, K2 = 4, K3 = 288, and 
K4 = 82944. In general one has the following result. 

Lemma 2 For any k > 1 one has 

2 fe - 1 fc!(fc-l)!(3£:-3)! 



{2k- 1)! 



Proof: We find a quadratic induction formula for the Kk , namely, for any k > 3 
one has 

(fc-2K = g( k ~ l ) ( ^ + ! ) ( k k _i 

Indeed, the left-hand side counts Wick contraction schemes together with the 
choice of an internal edge of the corresponding trees. If one cuts that dis- 
tinguished edge, the tree falls apart into two trees Xi, T<i. The numbering 
is unambiguous, if one decides that the cut edge goes from T2 to T\. Let i, 
1 < i < k — 2, be the number of internal vertices of T\. Choosing them ac- 
counts for the first binomial coefficient. The second binomial coefficient is the 
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number of ways one can choose the i + 1 external ip vertices of T\ among the 
inital k = ^. The third coefficient is for the choices of the k — i external tp* 
vertices in Ti- Note that the cut edge introduces an extra distinguished tp* leaf 
for Ti and an extra tp leaf for Now the given formula for Kk satisfies the 
quadratic induction because of the identity 

^ (3j)!(3fc-3j-3)! _ (3fc-2)! 

iK2i + l)!(2fc - 2j - l)\(k - j - 1)! k\(2k - 1)! U 

which holds for any k > 1, and which is an easy consequence of [25, Eq. 5.62]. 



Remark 2 One can also prove (6) using Clausen's 4F3 hypergeometric sum- 
mation formula. The k& 's are related to the Fuss- Catalan numbers of order 3 
(see [25, p. 347]). 

Now the rough estimate for the ^-clustering bound which comes from the 
aforementioned analysis of the tree level contribution, for even n > 2, is 

c n (A) ~ Cs^At- 1 ^ ~ Cs^At-^! 

by Stirling's formula and Lemma 2. 

It is therefore natural to make the following conjecture. 

Conjecture 1 The exists a constant c(J) > 0, such that for A > small 
enough, for any even integer n > 2, for any finite volume A G Z d , and for any 
assignment of the ft symbols, one has 



-— 1 1 

• n! . 



£ |(^ 1 (0),^ 2 (x 2 ),...,V tt "(x„))I, A | <c(J)"Af 
x2,...,x„eA 

Although we do not see a fundamental reason against this conjecture, due 
to technical difficulties inherent to the cluster expansion method we follow in 
this paper, we have been unable to prove so much. Related bounds of the n\ 
type have been obtained in the literature [20, 41, 13], but without extracting 
the optimal power of A. Using the method of this article, it is possible to 
extract this power, but at the cost of a higher power of n\. The difficulty is in 
obtaining optimal bounds which are uniform in n. Theorem 4 is a weakening 
of this conjecture. 



2.3 Free propagator decay 

Picking up the thread from §2.1, we let denote the off-diagonal part of J, 
so that J = J '(0)1 + J^. The matrix is also Hermitian and, because of the 
previous inequalities, has its operator norm bounded by ||J^|| < < J(0). 
Thus the Neumann series 
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converges, and provides the following random path representation for the free 
propagator 

1 1 (— l) p 

c(x ' y) = 7(0) 5x ' y ~ 7W t{ ^ y} J(x " y) + § 7 ^ TT 

x ^ J(x - zi) J(zi - z 2 ) ■ ■ ■ J (z p _2 - z p _i) J(z p _i - y) (7) 

zi,...,z p _i 

where the last sum is over sequences of sites in A such that, zi / x, z p _i / y, 
and Zj / Zj_i for 2 < z < p — 1. We assumed the function J : Z d — > C to be 
compactly supported. Let therefore ro > be such that J(x) = for all x £ Z d 
satisfying |x| > tq. We now have the following elementary lemma. 

Lemma 3 For any A C 7j d , and for any x, y G A one /ias i/ie uniform expo- 
nential decay bound 

|C(x,y)| <^oe" W)|x - y| 

w/iere 

J(0) 

k °-j,(j(o)-j,) >0 

and 

Hq = — log — — . 

r o \ J 

Proof: Given a nonempty A in Z, d , and x, y in A, let po = [ J X T , () y ^ J • First 
suppose that po > 1, so that |x — y| > r . Then obviously the first two terms 
on the right hand side of (7) vanish. Besides, if a term in the last sum over 
p > 2 is nonzero, then there exists a sequence of sites z l5 . . . ,z p -\ such that 
|x — zi | < ro, \z p -i — y| < ro, and |zj — Zj_i| < ro for 2 < i < p — 1. Thus 
l x — y| < P r o an d Po < Pi i- e -) P > Po + L Therefore at the level of matrix 
elements on has 



C(x,y) 



and bounding matrix elements by operator norms, one has 



from which the estimate easily follows, since po + 1 > ^ x ~ y ^ . Now if po = 0, one 
has |x — y| < ro and therefore using the full Neumann series and the coarse 
bound by the operator norm 

|0(x, y) | < ||C|| < 7M ^_ _ j^e-^ x ^1 
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and the last factor is obviously bounded by e w,r ° = ^51, so the lemma follows. 
■ 

This exponential decay trivially implies the I 1 property 

E^ |z| <^i(^) (8) 

for any jjl > and dimension d > 1, where Ki(d,fi) is some constant. 

2.4 The Brydges-Kennedy-Abdesselam-Rivasseau formula 

We recall in this section and the next one a list of well known combinatorial 
identities and inequalities and some basic results of cluster expansion that we 
will use to establish our bounds. 

The BKAR formula [15, 2] is a simpler and more symmetric version of the 
earlier Brydges-Battle-Federbush tree formula [14, 10, 12] in constructive field 
theory. This earlier formula itself is an improvement on the pioneering approach 
of Glimm-Jaffe-Spencer [23, 24]. 

Let us consider a finite set E ^ 0, and let us denote by E^> the set unordered 
pairs {a, b}, where a and b are any distinct elements in E. Of course \E^\ = 

' ' We will consider the space R E{2) of multiplets s = (si) leE (2) indexed 



2 

by pairs / £ E^ 2 \ and functions defined on a particular compact convex set 
ICe in this space. Let IIe denote the set of partitions of E. For any partition 
7r = {X\, . . . X q } in He we associate a vector v n = (^7r,/)i6B( 2 ) defined as 

v tt,1 = ^{3i,l<i<q,lcXi} ■ 

Now K.e is by definition the convex hull of the vectors v n , for ir £ He- It is 
easy to see that ICe affinely generates R s<2) . Indeed, let 6 be the partition 
entierly made of singletons, and for any pair / £ E^ let / denote the partition 
made of the two element set I and the singletons {a}, for a £ E\l. Then, the 
vectors vj — Vq, for I £ E^ form a basis of the vector space M. E 2) . As a result, 
the open domain Qe = ICe is nonempty, and ICe is equal to the closure £Ie- 
Let C k {CtE) denote the usual space of functions of class C k on the domain Q,e 
which, together with their derivatives up to order k, admit uniformly continuous 
extentions to the closure ICe = &e (see, e.g., [4]). 

Now a simple graph with vertex set E can be thought of as a subset of the 
complete graph E^ 2 \ A forest J is a graph with no circuits, and it is made 
of a vertex-disjoint collection of trees. Let $ be a forest, and let h = (hi)i G g 
be a vector of real parameters indexed by the edges I in the forest J. To such 
data we canonically associate a multiplet s($,h) = (s($,h)i) leE (2) in M. E 2> as 
follows. Let a and b be two distinct elements in E. If a and b belong to two 
distinct connected components of the forest J, then s(#, h)^ a h y = 0. Otherwise 
let, by definition, s($, h)^ a ^ = min h\ where I belongs to the unique path in 
the forest J joining a to b. We are now ready to state the BKAR formula. 
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Theorem 6 [15, 2] Let f G C\ E \~ 1 (n E ), and let 1 G R E(2) denote the multiplet 
with all entries equal to one. This is also the same as where 1 is the single 
block partition {E}. We then have 



/(1)= £ / 



where the sum is over all forests $ with vertex set E, the notation dh is for 
the Lebesgue measure on the set of parameters [0, lp, the partial derivatives of 
f are with respect to the entries indexed by the pairs belonging to and the 
evaluation of these derivatives is at the h dependent point s($,h). Such points 
belong to K-e- 

Note that the empty forest always occurs and its contribution is /(0) = 
/(^q). There are several proofs of this identity [15, 2, 16], but we believe the 
most natural and most easily generalizable is the one given in [3, §2]. This 
proof also most clearly shows the s($, h) belong to JCe- This point is important 
for the positivity of the interpolated covariance matrices in §4.1, and also when 
proving Lemma 5 via Lemma 4. 

We will now recall a lemma which, via the uniqueness of the Mobius inverse 
in the partition lattice, is a corollary of the BKAR forest formula. 

Lemma 4 Again let us consider a finite set E and let us denote by E^ the set 
of unordered pairs I = {a, b} in E. Let V( a ,6} be a collection of complex numbers 
indexed by E^ . Then 

1 tree 

Here g is summed over all simple graphs (i.e. subsets of E^ ) which connect 
E. We abreviate this property by the notation g ~» E. On the right-hand side 
the sum is on spanning trees T which connect E. The notation s(T, h) is as in 
Theorem 6. 

The following tree graph inequality, initially due to Brydges, Battle and 
Federbush (see [12, 10, 42]) is the basic tool we will need for the estimates in 
the first three regimes related to the 'high temperature' scenario. It is an easy 
consequence of Lemma 4. 

Lemma 5 Under the same hypotheses as in Lemma 4, let us assume that the 
numbers V\ satisfy, in addition, the following stability hypothesis: there are 
nonnegative numbers U a , for a G E, such that for any subset S C E one has 



aeS 
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Then the following inequality holds 



g~+E leg 



< e T,ae E Ua ^ Y[\Vl\ 



l e % 

X tree 



2.5 The cluster expansion for the polymer gas 

We now recall the basics of polymer gas cluster expansions. Any nonempty 
finite subset R C A is called a polymer. We denote by P(A) the set of all 
such polymers. We associate to each R £ P(A) a variable p(R) G C called 
the activity of the the polymer R. They can be collected in a vector p = 
(p(^))« e P(A)GC P ( A ). 

On the complex space C 13 ^ we consider the polynomial function Z defined 

P>0^ R lr ..,H p eP(A) k ^ 

for any p G C p ^. This function is usually called the grand canonical partition 
function of the polymer gas at finite volume A. It is well known (see e.g. [12, 17]) 
that the logarithm of Z\ can be written in terms of the following series 

log z(p) = £ 1 <f> T ( R i> ■ ■ ■ ' r p)p( r i) ■ ■ ■ p( r p) 

p>i P ' i?i,...,R P eP(A) 

with 

ff~>[p] 
ffCG 

where G is the graph with vertex set [p] = {1, . . . ,p} and edges corresponding 
to the pairs {i,j}, i / j, such that Ri n i?j / 0. The sum is over all spanning 
connecting subgraphs H which are identified with their edge sets. 

Note that by the so-called Whitney- Tutte-Fortuin-Kasteleyn representation 
(see, e.g., [53]), one can write the chromatic polynomial of G as 

P(G,x) = x c{H) (10) 

HcG 

where c{H) is the number of connected components of H. For nonnegative 
integer values of x, the quantity P(G, x) is by definition the number of proper 
vertex colorings of G with x colors. A good way to see the Ursell function 
(j) T (Ri, . . . , R p ) is as the coefficient of x in the chromatic polynomial P(G, x). 

The condition for the convergence of the series above is a well studied sub- 
ject. It can be expressed in terms of the following norm, depending on a pa- 
rameter a > 0, defined on the space C p ( A ) of polymer activities 

= su P Y \p(R)\ e am ■ 

xeA fteP(A) 
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The best result on this subject, essentially proven almost four decades ago 
by Gruber and Kunz [27] but largely forgotten, and then rediscovered very 
recently by Fernandez and Procacci [21] with a new proof, is the following 
theorem. 

Theorem 7 Let a > and let p denote an element ofC p ^ A \ Then the series 

= E i E ^ T ( fl i' ■ ■ ■ > r v)p(Ri) ■ ■ ■ p(R p ) 

p>i P ' fli,...,flpeP(A) 

is absolutely convergent in the closed ball \ \p\\ a < e a — 1. The function f is 
analytic on the open ball \ \p\\ a < e a — 1 and satisfies 

exp/(p) = Z 

for any p with \ \p\\ a < e° — 1. 

In fact, one can extract a more precise result (see [21, p. 132]), when p > 0, 
i.e., when the polymer activities p(R) are real and nonnegative. 

Theorem 8 If p > and ||p|| a < e a — 1, then for any Rq 6 P(A), we Ziaue the 
estimate 

E i^ T (^' ■ ■ ■ > ^)i-p(^i) • • • p( R P ) ^ ea|Ro1 - 1 • 

p>i P ' Ri,...,-R P eP(A) 

We will use this theorem when i?o is a singleton, namely, when Rq = {z} 
for some z G A. In this case one has the identity 

4> T (R ,Ri,...,R P ) = (-r). ( f> T (R l ,...,R p ) (11) 

where r is the number of indices q, 1 < q < p, such that z G Indeed, if G 
is the intersection graph on {0, 1, . . . ,p} defined by the collection of polymers 
Ro, Ri, . . . , R p , then the restriction of G to the set formed by the vertex and 
its neighbors is the complete graph on r + 1 elements. This property and the 
representation (10) gives a very easy proof of the reduction formula (11). A 
trivial consequence of (11) and Theorem 8 is the following lemma which will be 
used repeatedly in the sequel. 

Lemma 6 For nonnegative polymer activities p{R), R 6 P(A), such that \ \p\\ a < 
e a — 1 we have the bound 

su pEi E i{z e u^RJ 

zeA P>i 1 fli,...,flpeP(A) 

x\cf T (R 1 ,...,R p )\.p(R 1 )...p(R p ) <e a -l. (12) 

Remark 3 The important point to note here is that the bound is uniform in 
A G Z d . 
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Since we are not interested in optimal bounds we will choose hereafter a = 
log 2 and we will denote the norm ||p||io g 2 simply by ||p||. So in what follows 
we will use the norm 

Ml = sup \P(R)\ 2 m (13) 

xeA i?eP(A) 

with the condition ensuring absolute convergence of f(p) being 

< 1 • (14) 



3 The large mass, small interaction, and large self- 
interaction regimes 

3.1 The Mayer series representation for the truncated correla- 
tions 

Given a source specification (x,, (Ji)ig/, we will consider the perturbed partition 
function 

Z A (a) = [ Dil>*Dil> e -^A^A) TT(i + a ^( x .)) 
J ^ A tei 
and the truncated correlation function given by 

B\i\ 

<(^(xi)) ie/ )J = Log Z A (a 



a=0 



As mentioned before, Log Z A {a) is analytic in a small polydisc D A around 
a = 0. 

Now let us introduce the normalized single site measure 

du(z*,z) = -le^WN 2 -!! 2 ! 4 ^ d%z 

on C, where 

N= f e -A0M 2 -$\*\ 4 Mz d%z . 
Jc 

Clearly Z A (a) = N^Z A {a) where 

Zk{oc) = / II cfc/(V*(x),V(x)) e -^{x, y}e A( 2 ) WW TJ(i + a .^( Xi )) 
•^ cA xeA is/ 

with 

WW = J(x - y)V*(x)V(y) + J(y - x)y,*(y)y,(x) . 

bmce Z A (a) = AA |A| Z A (a) with > 0, one has the analyticity in the polydisc 
D A of Log Z A (a) = — |A| In N + Log Z A (a) as well as the identity 

MH*i))iei)I = n ^ ' Log Z A (a) 
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On the other hand, we can also write 



e -£ {x , y}eA (2)J { x,y}(iW = Yl [l + (e-WW _ 1) 

{x,y}eA( 2 ) 

= E II (e- J ^>W-l) 

gCA( 2 ) {x,y}eg 

where the sum is over all simple graphs g on the vertex set A. Using this equa- 
tion and also expanding the product of the (1 + ajV^( x «)) one easily obtains, 
after reorganization according to the connected components of the graph g, the 
following polymer representation for Z\(a). Namely, one has 

**<«)- £3 E »{*^ }**>,»).. .^a) (15) 
p>o 1 fli,...,flpeP(A) k ^ 

where the polymer activity of a polymer i? is defined as 

I HdvirW^Wil^fr)) H (e"W*)-l) (16) 
JcR xeR ieJ {x,y}e g 

with the notation ^ = {16 J|xj 6 i?}. Let p(o:) denote the activity specifica- 
tion (p(.R, a))^ 6 p(A). Using the notation and definitions of §2.5 one obviously 
has Z\{a) = Z(p(cx)). 

Now suppose the condition ||p(0)|| < 1 holds. Then on a small polydisc 
D' k the hypothesis ||p(a)|| < 1 will also hold. On the polydisc D\ n D' A both 
the functions f(p(a)) and Log Za(cc) are analytic and exponentiate to Z\{a). 
By connectedness the difference f(p(a)) — Log Z\{a) is a constant in 2iirZ. 
However this difference takes a real value at ex = 0, and therefore it must vanish. 
Indeed, the hypothesis on the J(x) function implies that // Xiy }(Y>) is real and 
therefore, by (16), the activities p(R,0) are also real. Finally, by the definition 
in Theorem 7 it follows that /(p(0)) £ I. As a result of these considerations, 
on a small polydisc around a = one has the equality 



/(p(a)) = Log Z A (a) 



thus 



ct=0 



= E^ E m,-.^) 
Hi 



p! 

p>o fli,...,flpeP(A) 



■p(i?i,a) . ..p(R p ,a) 



a=0 
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= Ei E <p t (Ru-..,r p ) E 1 



p>o P ~ fli,...,flpeP(A) 



h,...,I p Cl 



Iq disjoint 

uj, = / 



n 



0^9 I 



■p(R q ,a] 



a=0. 



Comparing with (16) it is apparent that the effect of the a derivatives is to force 
the summed over J to be equal to I q . Hence, as an immediate consequence of 
Theorem 7, we have the following statement. 

Proposition 1 Define for any subset J of the source label set I, and any poly- 
mer R £ P(A), 



' Xj G R J lorJ/0 



E 



Provided the source-free condition 



{x,y}eg 



sup 



xeA 



KeP(A) 



holds, one has the absolutely convergent series representation for all truncated 
correlation functions 

(W»'(x,)k,)L = 

r>>l P- H, flp£P(A) k 9 J _ =1 



/i,...,/ p c/ 



Note that with this new definition, polymer activities p(R, J) do not have to 
contain all the sources localized in sites belonging to R. 



3.2 Estimates for a single polymer activity 

Consider a polymer activity p(R, J) as defined in Proposition 1. By moving the 
sum over connecting graphs inside the integral, one has the estimate 

(ni^(*)i)* e n (e- 7 <-> w -i) 

Vie J / |g^fi{x,y}eg 

Now for any subset S C R and for any fixed field tp, by the argument given in 
§2.1, one has 



E J {x,y}(V0 
{x,y}es( 2 ) 



<J^^|^(x)| 2 . 

x<=S 
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As result, Lemma 5 implies, for any given field configuration ip, the inequality 



< 



II (e-^-> W -l) 

g-+-R{x,y}eg 

e J ^ R \^)\ 2 II (2|J(x-y)||^(x)|^(y)|) . 



| T h { x, y}e S 



Therefore 

\p(R,J)\ < i 



Vi G J 

X,; G 



or 



} 2lfi| - 1 e n 

' J rr . (J f 1 ,- ,-7- 



{x,y}eX 

X tree 



X I1 / ^(fWi(x))e^ Wx)P Wx)r (x) 
xeir " 



C 



(17) 



with m(x) = cj(x) + d<j:(x) where cj(x) is number of source labels i G J such 
that Xj = x, while %(x) is the degree of the vertex x in the tree 1. The 
integrals in (17) are now estimated thanks to the following lemma. 

Lemma 7 For any integer m > 0, and under the hypotheses of the Introduc- 
tion, we have 



J dv(z\z)e J ^\z\ m < 2 ^J(O) + (J(0) 



m + 2 1 

- J # ) — ^m!2 (18) 



as well as 



/ dv(z*,z) 
Jc 



(19) 



Proof: We have 



du(z*,z)e J ^ 2 \z\ m = -1 / e -( J (°)-^)kl 2 -|l^l 4 | z |™ d^z dQz . 
In order to obtain a lower bound on the denominator ftf, we write 

oo oo 

N = J c e- J WW 2 --^dMz d%z = 2vr I e - J W-^rdr = tt J e'^'^dt 



oo 

f -MVi 
= 7re a e 



t+J(o) 



[2 J(Q) 2 f 

dt = 7r \lx e A J 



e 2 dq 



J(O) / 2 



Now recall Birnbaum's inequality [11] for Mill's ratio, i.e., essentially the 
erfc function: 



/ 



e 2 e 2 dq> 



x 
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for x > 0. One can simplify this to 



oo 

sL f _ai , 2 
e 2 / e 2 dq> ^= > 



x 



+ \/4 + x 2 x + \/4 + Ax + x 2 x + l 



The latter applied to x = J(0)y | readily provides the needed bound 

M > • (20) 

J(0) + ^\ 

Now for the upper bound on the numerator, we have two possibilities corre- 
sponding to the two estimates in the lemma. 
1st option the bound using the Gaussian part: We write 

I e -(J(°W*)l*l 2 -$l*l 4 \ z \m Mz d Q z < f e -{m-W \ z \m Mz d c Sz 

./; Jc 



= 27T j e -(m-^y rm+ 2± = 27r(J(0) _ j^-^r . 



We therefore only need to show the elementary inequality < 1 in 

order to complete the proof of (18). This can be done easily by induction using 
the well-known properties of Euler's gamma function (see e.g. [5]). The quantity 
-i=r (™±2) i s equal to 1 and ^ < 1 for m = 0, and m = 1 respectively. 
Besides, when going from mtom + 2, this quantity changes by a factor of 



1 m + 2 1 I 1 

x — — = -«/! + r < 1 



• s /(m + l)(m + 2) 2 2V m + 1 
and the desired inequality propagates. 

2nd option the bound using the quartic self-interaction: We write 



f e ~(J(0)-J^)\z\ 2 -±\z\* | z |m ^ < /" e -| k| 4 | z |m ^ ^ 

./:• 7: 

00 , m + 2 



We are left with showing T < 2m! 3 for any nonnegative integer m. By 

Dirichlet's multidimensional extension of the formula for the beta integral [5, 
Theorem 1.8.6], we have 



/m + 2\ = ( m + 1 y[ (tititsU) ^ dt 1 dt 2 dt 3 dt 4 

V 4 / Jt t >o,Y,u=i 

to — 2 

< (m + 1)! / ( — J * dhdt 2 dt 3 dt 4 

7t i >o,E*i=i V 4 / 
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by the arithmetic vs. geometric mean inequality. Thus 

r (^j^j 4 < (m + 1)! x 1 x 4~ m + 2 < 2 4 m! 
for m > 0. Indeed, 

4 m = (1 + 3) m > 1 + 3m > ( m +^ . 

■ 

From the raw estimate (17), the previous Lemma, the inequality 
m(x)! < 2 Cj ( x ) +d ' I ( x ) Cj (x)!%(x)! 
and the trivial relations 

£ Cj (x) = |J| , £<fc(x) = 2|i2|-2 

xS-R xS-R 

one easily derives the following basic bounds on single polymer activities. 

Lemma 8 Gaussian estimate: 

One has the bound 

\ R \ 

x(j(o)-j^)- 2 i r i-^ +i x n^w^x e n^w 1 n i j (*-y)i- 

xeR X~*R xeR {x,y}eT 

X tree 

Lemma 9 Quartic estimate (a.k.a. domination): 

One has the bound 

xA*t+Ix[] C j(x)^x ^ n%(x)!i [] l«/(x-y)|. 

xei? s^fl xei? {x,y}el 

% tree 

3.3 The large mass and small interaction regimes 

In this section we will use the estimate of Lemma 8 which allows us to treat at 
once both the cases when J(0) is large or when is small. Our first task is to 
check the source-free (i.e. J = 0) condition for the applicability of Proposition 
1. Let z be a site in A. We need to bound 

A = H^R} \p(RM 2 |R| 

i?eP(A) 

< £ l{z € R}1{\R\ > 2}2 4 l*l- 2 (j(0) + Jf] 



x 



(j(o)-j^i +1 x^[]^ n l J ( x -y)l 



s~>« xei? {x,y}ex 

T tree 
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We now condition this sum on the size m = \R\ > 2 of the polymer, and 
introduce a sum over labelings of the sites in R by the fixed set of indices [m]. 
Namely we introduce in the sum the identity 

1=-—— - l{xi =z}l{ Xi distinct} l{i? = { Xl ,...,x m }} . (21) 

xi,...,x m eA 

The next step is to use this artifice to transport the summation over trees T 
on the variable set R to a sum over trees t on the fixed set [m]. Then, partly 
releasing the second condition in (21), and eliminating R we obtain 

^E^ E «*. -.>*"-' + 

m>2 v ' xi,...,x m eA \ / 

x(j( )-^)- 2m+i x n^) !l n (%^ 3 }|j(x l -x 3 )D 

t tree 

with the obvious notation c#t(i) for the degree of z € [to] in the tree t. Now the 
sum over the locations Xj , starting with the leafs and then progressing towards 
the root 1 G [m], is easy and gives a factor J^~ l . The sum over the tree is 
done using the following easy lemma. 

Lemma 10 We have the bound 

£ II d i {i)\<2^\m-2)\ . 

t~>[m] i£[m] 
t tree 

Proof: We have by Cayley's Theorem (see, e.g, [54, Theorem 5.3.4]) which 
counts labelled spanning trees with fixed vertex degrees 

Efi*w - E Y. }W 

t~»[m] i=l il,...,d m >l t— [m] <* i=l 

t tree Ed 4 =2m-2 t tree 

= „ (m-2)! ft 

IE.(*-i)iii*' 

Ed i= 2m-2 

= (m - 2)! £ d 1 ...d m . 

d\ , .. .,d m >1 
Edj=2m-2 

Using the arithmetic versus geometric mean inequality 

"2m -2" 



di...d m < 
as well as 



I m 

< 2 m 



m 



E i=(^_-, s )s2 



2m-3 



Ed i =2m-2 



21 



the lemma follows. I 
Thanks to the coarse bound dt(i)! 5 < dt(i)\ and the last lemma, we now 
have 



A < £ 2 7m ~ 5 J™" 1 (j(0) + (J(0) - 

m>2 V / 



2m+l 



Therefore, as soon as the condition 



2%. (j(0) + ^ x 



(J(0) - J # )* 



holds, one will have 



(22) 



l|p(-,0)ll<A< 



2^.(j(0) + y|) 
( J(0) - J # ) 3 



(23) 



Clearly, if we take either one of the limits J(0) — > oo or — > 0, the condition 
(22) will hold and ||/o(-,0)|| can be made arbitrarily small. Note that crucial 
to this last fact is m > 2, i.e., the absence of single-site polymers, also called 
monomers. 

We are now in a position to tackle the ^-clustering estimate (4), where 
the source label set is / = [n\. Assuming either J(0) is large enough or is 
small enough to garantee conditions (22) and ||p(-,0)|| < 1 hold, we can use 
Proposition 1 as well as Lemma 8 to write 



£ KV ai (o),V' fe (x2),...,^(x„))] 

x 2 ,...,x n eA 



xi,...,x n eA 



m 

p>l 1 flj flpSP(A) 

,/pC[n] 



n 



q=l 



Vt G J, 

X,; e i?„ 



1^1 > 2 
or |J,| > 1 



x2 3\R q \ + M-2 



n^(x)i^x ^ n°M x ) !l n i j ( x -y) 



xei?„ 



Tq tree 



{x,y}eT„ 



(24) 



The first step is to push ll{xi = 0} through the sums over p, the i? g 's and the 
iq's. We then bound it by the coarser condition ]1{0 € Ug =1 i?q}. The second 
step is to push the sums over the Xj's inside the appropriate (i.e., as dictated 
by the choice of the J ? 's) bracket factor. The sums over the source localizations 
Xj are therefore bounded with the knowledge of which polymer they belong to. 
This rests on the following lemma. 
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Lemma 11 Using the notation c y (x) = \{j, 1 < j < k\yj = x}|, /or any 
polymer R, and for any power (5, we have 

II c y( x ) ! " ^ 2|i?l+fe_1 x &! max{/3,1) • 

Proof: Summing first over multiindices c = c(x) x£i r and then over the se- 
quences (yi, . . . , yfc) for which c y = c, we have 

x ri c y( x ) ! "= fc! e n^)^ 1 

yi,...,y fc eR x£R c,|c|=fc xSi? 

where we denoted by |c| the length X] x eR c ( x ) °^ * ne m ultiindex c. We then 
bound the product by /c!^ -1 if (3 > 1 and by 1 otherwise. We also use the trivial 
bound 



^ 1= | \ R \~^ 1 + k ^j < 2 \ R \+ k - 1 



c,\c\=k 

and the result follows. 

As a result of this lemma we have for every q, 1 < q < p, 



e A *| H >< n c ^w !l ^ 2l ^ l+|/9hi x \w • w 

(xi) i€J ,eA J » * q xei? 9 

Calling 2B the sum to be estimated on the left-hand side of (24), and using the 
relation J2q=i \Iq\ = n in order to pull out some factors from the sums, we now 
have 



2H < 2^.(J(0)-^)-f^i £ WeU^M 

P>1 Ri,...,flpeP(A) 

*\f{R u .^R p )\l{ Jxl^lIx-.-xl^l! 



p 

«n 

9=1 



11 



\Rq\ > 2 

or |J g | > 



\R q \ 



:(j(o) - j^-w 1 x ^ n^w !l n i j ( x -y) 



1 (7 tree 



We now insert the decomposition 

A i R n^\ ] = H\R q \>2} + A J^i > (2«) 



or |7 9 | > 1 J g [ and |I g | > 1 

so each bracket factor takes the form A q +B q where the nonnegative numbers A q 
and B q correspond to the first and second conditions of (26) respectively. Note 
that since the I g 's form a disjoint decomposition with possibly empty subsets 
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of the set [n], the number of g's for which B q / is bounded by n. Hence, for 
any number 7 such that < 7 < 1, we can write 



f[(A q + B q ) = f[(A q + B q )x [] A q 

q=l 9=1 9=1 

B q ^0 B q =0 

< f\{r 1 A q + B q )x fl A q 

9 = 1 q = l 

Bq^O S 9 =0 



g=l 



9=1 



9=1 

where B g means we now forget about the condition \I q \ > 1. Leaving the 
appropriate choice of 7 for later, we now have 

2U < 7- n 2¥.(J( )-^)-i^i Yl HOelf q=1 R q } 

p>l P' fi lr .,B p £P(A) 
I 1: ...JpC[n] 

x|^(fi 1 ,...,R„)|j{ ^T^' }x|'i|U...x|/ p |! 



9=1 



x 



(j( ) - x ^ n^w 1 n i j ( x -y)i 



To tree 



{x.yjeX, 



+jt{\R q \ = 1} x 2 ( J(0) + ^ j x (J(0) - J^)- 1 



We now perform the sum over the I q s 



E 1 

h,...,I p C[n] 



I q disjoint 
Ul q = [n] 



E 



k 1 ,...,k p >0 
fclH f-fcp=n 



n 

/ci • • • /c 



x |/i|! x ••• x \I p \\ 
^ x k\\ . . . kp\ 



n +p — 1 
p — 1 

< n!2 n+p_1 . 



(27) 
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Thus 

p>i P ' fli,...,flpeP(A) 
H{0€U^ =1 i2,}|^ T (12i,...,12 p )|.e(i2i)...e(i2 p ) (28) 

where the nonnegative polymer activities £>(•) are defined by 

J(o) + ■ F\ 

q{R) = 4j.l{\R\ = 1} x 



J(O) - 

-2|fl|+l 



> 2} 2 4 l*l~ 2 (j(O) + Jfj (J(O) - J^)~ 



n^w 11 n i j ( x -y)i 



The norm is estimated in the same manner as we did for ||p(-,0)|| at 
the beginning of this section. Namely, we find 

J(o) + y| 

1 |p| 1 < 87 x 



J(0) - J # 

+ 2 8m - 5 j™- 1 ( j(o) + yf) (j(o) - j^y 2m+1 



m>2 



By taking J(0) large or small we will ensure that 

2 8 ^. (j(0) + yf^j 1 

(.7(0) - J^) 2 -2 ■ 
and also, summing the geometric series, that 



7(0) + ^ 2*J+(j(0) + J$) 



llel|£87X 7(oF^ + (.7(0) -j,). (29) 

holds. 

1st case: When J(0) becomes large, the right-hand side of (29) approaches 
87. So we simply choose 7 = ^ < 1 and Lemma 6 together with the previous 
estimate (28) complete the proof of Theorem 1. 

2nd case: When becomes small, the right-hand side of (29) approaches 



87 ^1 + ■ Therefore, we choose 



7 = t r- < 1 



!6 + ^ 



and now Theorem 2 follows. 
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3.4 The large self-interaction regime 



This essentially is an encore presentation of Section 3.3. The difference is that 
we now use the estimate of Lemma 9. Following the same line of argument as 
in §3.3, we therefore successively have 



A = Yl H^R} \p(RM 2 |R| 



ReP(A) 



< M» e R}1{\R\ > 2}2^-f (j(0) + Jfi 

x Y n^( x ) ! " n i j ( x -y)i 

then, introducing the labeling and using Lemma 10, 



\ \R\ 



A<Y 2 j # m "^.A~ m+ 5.J^- 1 ( J(0) + j| 

m>2 V y 



For A large we will have 



and 



So we now have, 



2^.(J(0) + Vt / , 
- 2 



||p(-,0)||<2"j # A-f (j(0) + ^) <1 



w = Y i ^ 2 ( X2 )' • • • > ^ n ( x «))aI 

x2,...,x„eA 

< E ^^E^ E 



xi,...,x n eA 



p>l ^ Hi,.. .,ftpgP(A) 

'i 'pCW 



.:/■, D n> m ii J disjoint 
Ul q = [n] 



\Ri,...,Rp)\l 



n 

? =1 L 



1^1 > 2 
or |J„| > 1 



x2 ¥i^i+fi^i-f ( j( ) ■ /A 



n c/9 ( X )!ix ^ n^( x ) !i n i j ( x -y)i 

{x,y}eT„ 



'X o tree 



(30) 



Following the same steps as before, including the introduction of 7 6 (0, 1], we 
arrive at 



2H < nl.T-^^.A-fx^i ^ 



P>1 Ri,...,R p 6P(A) 



26 



1{0 G U p q=1 R q }\<p T (R u . . .,R p )\.<;(R 1 ) . . . s(R p ) 
with the nonnegative polymer activities 



,(R) = 2 4 7 .1{| J R| = 1} x A-H J(0) + A /^ 



\R\ 



:X -\r\h y: II i^( x -y)i- 



57 K xeR {x,y}6l 

X tree 



Their norm is similarly bounded by 

Ml < 2§ 7 x ^1 + J(0)^ 

+ J] 2% m - 1 ij™- 1 \~ m+1 2 (j(0) + W^j 

m>2 V / 

For A large enough we will have 

2f^(l + J(0)y|) i 
A - 2 

and therefore also 

— \ / i — \ 2 

2\ 33 , . 3 / , . . /A 



Ml < 2n x ^i + j(o)y - j + 2-j^a-5 ^j(o) + y - 

9 

The latter expression approaches 227 as A becomes large. Therefore, choosing 
7 = -4r and applying Lemma 6 completes the proof of Theorem 3. 

2T 

4 The small self-interaction or near-Gaussian regime 

4.1 The cluster and Mayer expansions for the truncated corre- 
lation functions 

We will follow a line of argument similar to §3.1 in order to obtain a convergent 
series representation of the truncated correlations which is adapted to the small 
A regime to be considered in the remainder of this article. Using the notation 
of §2.1 and §3.1 we have Z\{a) = (det J) -1 Za(o:) where 

Z A (a)= [ ^ c (f,^) e-i^Al^)l 4 []fl + a ^(x 8 )) . 

Since J is positive definite, we again have that Log Z^{pt) is analytic in a small 
polydisc Da around a = and 

(tyH*i))iei)A = n ^ ' Log Za(<*) 
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The next step is to write an expansion for Z\(ct) similar to (15). This step is 
usually called the cluster expansion in the constructive quantum field theory 
literature and was introduced by Glimm, Jaffe and Spencer [23, 24]. These 
expansions which have been simplified and improved over the years by a small 
group of experts, are not so well known in the wider mathematical community. 
In what follows we will try to explain the method in detail, on the simple case 
of the lattice field theory considered in this article. We also use one of the 
more recent technical implementations based on the BKAR forest interpolation 
formula of §2.4. This is a kind of combinatorial Taylor expansion with integral 
reminder which interpolates between a 'complex' fully coupled situation and a 
'simpler' fully decoupled one. In the present case this decoupling expansion will 
be applied to the Gaussian measure, since it is the only feature preventing the 
random variables of different lattice sites from being independent. 

Before introducing the decoupling expansion for the Gaussian measure, we 
need a preliminary expansion of the self-interaction term e~ * ^ x e* l^( x )l . This 
is a matter of convenience and spares us the division by the amplitude of triv- 
ial polymers when deriving the polymer gas representation. This is especially 
useful in view of the forthcoming derivatives with respect to the coupling con- 
stant A which one would rather have act on products instead of ratios. This 
preliminary expansion consists in writing 

e -| E xeA W*)l 4 = J] [l + (V^MI 4 - l)' 
xeA 

= e n («-*'* x)|4 - 1) • 

TcAxeT 

Then for each x G T we write 

l 

e-^WI 4 _ 1 = J dt x (-^(x)| 4 ) e~M^)\\ 
o 

As a result, one has 

I ^ I 

Z A (a) = (-7) / £ I d^cir^) 

V 7 [0jl] T J C A 

^n^ + ^'fc))) (il I^WI 4 ) e-iExa'*WI 4 (31) 

where dt denotes the Lebesgue measure 11x6 r ^*x- 

We are now ready to apply Theorem 6, using E = A, as follows. For any 
multiplet s = (s/) ;gA (2) in the closed convex set /Ca, we replace the covariance 
C by a modified covariance C[s] defined by 

C[s](x,x) = C(x, x) for any x G A, 

and 

CM( x >y) = s {x,y}C(x,y) for x / y in A. 



28 



Clearly this new matrix is also Hermitian. Moreover, the following key positivity 
property follows from the previous definitions. 

Lemma 12 If s £ K\ then C[s] is positive definite. 

Proof: For s G /Ca one can find an expression s = J2j=i w j v ^j with the Wj's 
nonnegative and summing up to 1, and the TTj some suitable partitions of A. 
Using notations similar to §2.1, for any vector ip G C A we have 

<V,W> = £ V>*(x)C[ S ](x,y)V>(y) 

x,yGA 

= J>*(x)C(x,x)V(x) 

xeA 

k 



+ E^' E V'*(x)K J ){x,y}C(x,y)V'(y) 



j = l x.yeA 



k 



= E 



w 



^V*(x)C(x,x)^(x) 



.xeA 



x,y£A \Xewj 

j = l XETTj 

where ipx( x ) = 1{ X G X}tp(x.). Therefore the result is nonnegative since C is 
positive definite. Now suppose < tp,C[s]ip > vanishes. Since Ylj=i w j = 1> one 
can choose jo such that Wj > 0, and then have for every block X G irj that 
< ip x , CVx >= 0, i.e., that ip x = 0. Thus V = Ex 67% i>x = 0. ■ 
As a result of this lemma the corresponding Gaussian measures d/U^[ s ] are 
well defined. The last Gaussian integral in (31) therefore becomes a function 

ns)= j cA dnci.](r,ii>) ^n^+^fc))) (^ni^ x )i 4 ) e ^ E ^ txWx)|4 

of s G /Ca to which one can apply Theorem 6. The well known rule for com- 
puting derivatives of Gaussian integrals with respect to the covariance matrix, 
namely as a Laplace type operator acting on the integrand (see e.g. [22, §9.2]), 
immediately implies the following representation: 

- e B) 1X1 L * £. L "V- <w-*> 

on A 

,/e5 / \iei J VxeT / 
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Here we denoted by A; the operator of Laplace type given, for any unordered 
pair of distinct sites x and y by 

A {x , y} = ^(*,y)^y^y + C(y>*)a^) (32) 

using the standard J| and J= vector fields of multivariate complex analysis. 
The differential operators act on everything to their right. 

Now, let 7r be the partition of A into connected components of the forest 
J. An important feature of the definition of s($, h) is that it vanishes between 
components. This implies the componentwise factorization of the Gaussian 
integral. One can also factorize the different combinatorial sums involved, i.e., 
those over the sets T# = TDR, as well as the ones over the trees Tr connecting 
each R £ ir and which together make up the forest In sum, one has 

Z A (a) = n^' a ' A ) ( 33 ) 
7ren A Re-* 

where for each polymer R we defined the polymer activity 

cb(fl,a,A) = E E (-j) I d * [ dl 

T tree 

vex I \ieJ I VxeT / 

(34) 

again with Ir = {i £ /|xj G R}. Note that the covariance C which is used here 
is the restriction to R of the one defined on A as the inverse of J. It therefore 
retains a slight dependence on the volume A which contains the polymer R. We 
nevertheless suppressed it in the notation for better readability. 
We now introduce in (34) the decomposition 





' \R\ 


= 1 and 




' \R\ > 2 or ^ 


1 = 1< 


T 


= and 


M 


T > 1 or 






J = 




{ \J\>1 J 



and break Co(-^> a , A) accordingly as a sum of two contributions. It is easy to see, 
because all sums and integrals become trivial and also because even for a single 
site the Gaussian measures are normalized, that the first contribution reduces to 
= 1}. Therefore, by only keeping track of polymers for which the second 
contribution is selected, one can rewrite (33) as a polymer gas representation 
similar to (15), namely 

*a(«) = E h E A 1 dlfoint 6 } C(*i, «, A) • • • «Rp, «, A) (35) 

P>0 F ' J Ri,...,i? p eP(A) k J 
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with the polymer activities defined as 

' \R\>2 or 
|T| > 1 or 
|J|>1 



c(R,c*,\) = j2 E E 1 

rcRJci R ^ 

\teJ J J ^ r J ^ JcR 



4 



(u a ) (u^w) fn i^wi 4 ) "■^^ • (36) 

\iel J \i<=J J VxeT / 

Now the same line of argument leading up to Proposition 1 shows the following. 

Proposition 2 Define for any subset J of the source label set I, and any poly- 
mer R £ P(A), 



™^^} T ¥«H) IT 'l 

% tree 

{wA fn^wl fn i^wi 4 ) e-^^Mr (37) 

V/eX / VieJ / VxeT / 



' > 2 or 
T| > 1 or 

1^1 >1 



Provided the source-free condition 

sup l{x G i?} |C(^,0,A)| 2l R l < 1 (38) 

xeA i?eP(A) 

holds, one has the absolutely convergent series representation for all truncated 
correlation functions 

<(^)W)I = E^ E 4> T (Ri,-..,Rp) 

p>l V' fli,..,Jip£P(A) 
I 1: ...,IpCI 

x4 J ^r 7 nt )nc> ? ,/ ? ,A). (39) 

9 J q=l 

Remark 4 Although this is not really necessary for the proof of Proposition 
2, note that the source-free activities C(-R, 0, A) are real-valued. This is because 
of the hypothesis J(— x) = J(x)* which implies that the differential operators 
A{ X y} preserve the real-valuedness of functions. 

Remark 5 In the constructive quantum field theory literature, an equation such 
as (35) would be called a cluster expansion for the partition function Z(a) which 
is written as a sum over collections of disjoint polymers. An equation such as 
(39) involving sums over collections of polymers with the coefficient T would 
be called a Mayer expansion for the truncated correlation functions. 
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We now go back to the setting of Theorem 4 and also restore the A depen- 
dence in the notation for correlation functions. In order to extract the factor 
X N in the clustering estimate, we will perform an additional Taylor expansion 
of the connected correlation function. We write 

\k f a \ k 

<^(X 1 ),...,^-(X B )>3[ IA = £ ¥ - (^(xi),...,^"(Xn))lJ =o 


However, by Lemma 1 and the hypothesis n > 2(N + 1), the Taylor polynomial 
vanishes and the right-hand side reduces to the integral remainder, i.e., 

(^(xi),...,^(Xn))I A = 

/ du lW~^r (^) iV ^ 1 ( Xl )'---'^"( x -))I^ ■ 



So far we only used the fact that correlations as C°° functions of A on [0, +oo). 
However, the next step is to use the representation (39), and also to differentiate 
it, term by term, N times. This will require some estimates, in order to justify 
the convergence criterion (38) as well as the following outcome of term by term 
differentiation 

W,»'(x 1 ),...,«'"(x„))I, A = 

/ £ * T w....,« P ) 

i V >' p>l 1 -R 1 ,...,H P SP(A) 

/l,...,/pC[n] 

for n > 2(N + 1). Finally we will use this last representation as input for the 
^-clustering property. 

4.2 The estimates 

This section will provide the necessary estimates for the proof of Theorem 4. 
Many of the ideas used in these estimates originated in the work Glimm-Jaffe- 
Spencer (see e.g. [22, Ch. 18], [44, §111.1], and [3]). We will first suppose that 
the J function governing the Gaussian measure is so chosen that the constant 
Kq appearing in Lemma 3 is equal to 1. Later in §4.2.4 we will get rid of this 
restriction by a simple scaling transformation on the field variable tp. 
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4.2.1 The bound on a single polymer activity and its derivatives 

The basic quantity we need to bound is {-§^) M C(R, I, u\), for any integer M > 
0, polymer R, index set I € [n], and u € [0, 1]. From (37) one easily gets 

j \ M ( w . T ( \R\ > 2 or 

£) «^) = »(^h)E» mil 



or 

TC-R I 171 > 1 



du] v ' ' ; [ Xj G 

ITI I / \ n ' - 

L. j 

k=0 



E E / d * I dlUty*) 1 ^ 



X tree 

where X refers to the remaining Gaussian integral, namely, 



vzel / Vie/ 



A; 



^11 I^WI'J [U^^j e-*^rW)\* . (42) 

When one performs the derivatives coming from the A; operators, this further 
splits into 

^=E J p 
p 

where, avoiding excessive formalization, we denoted by p any derivation proce- 
dure including the detailed information as to which specific tp or ip* factor has 
been destroyed, and by which or -^p- operator. We will use a bound of the 
form 



\1\ < (max\Ip\jJ2l 



Note that a term J p has the form 



Ip — Cp.Lp 



J] {^(x) m Wf(x) m *W}| e-TLa *xW*)l 4 



where the m(x) and m*(x) are some local multiplicities and L p is a product 
of factors — A«£*x produced by the derivatives which acted on the exponential. 
Finally, Cp is a product of propagators corresponding to the edges in the tree 
T. It depends on whether, for each A/, one choses the C(x, y) or the C(y,x) 
term. In any case, Lemma 3 together with the assumption Kq = 1 implies 

\Cp\ < e -«£ {x , y}es l*-yl . 
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We will use the notation \m\ = S x eR m ( x ) ^ or ^ ne total multiplicity of a mul- 
tiindex such as (m(x)) xe #. Using the positivity of the interaction, and in 
particular the positivity of A, one has 

< e-xExa'xI^WI 4 < i . 
Therefore, by the Cauchy- Schwartz inequality 

d ^c[s { %h)p*^) ni^(x) m(x) v*(x) m * (x) } e -TE xe T*x^(x)i 4 

xe_R 



i 




where m(x) = m(x) + m*(x). 

From now on we also impose the hypothesis < A < 4, which implies 
\Lp\ < 1. We now invoke the following classical lemma of constructive field 
theory [24, 20] also called the principle of local factorials, in order to bound the 
last Gaussian integral. 

Lemma 13 For any q, and collection of sites zi, . . . , z q and w±, . . . , w q in R, 

one has 

<KlY[n*(x)\ 

where n*(x) denotes the number of indices i, 1 < % < q, such that Wj = x and 
K^K^no) of (8). 

Proof: By the Isserlis-Wick Theorem we have 

j cR d ^C[ s (T,K)] C^*, VO V>( z i) • • • ^( z gW ( w i) ■ • • ^* ( w <?) 

< £ n|(7(z,,w CTW )| < ^ e -w»E?=il-i-w.( l 

(766,i=l ff£6, 

because the decoupling parameters in s(T, /i) are between and 1, and also 
because of Lemma 3. Now 

e -MoELil z i-w CT(i) l = e -wELiK-^l ^ l{Vi,w ff(i) = Ui } 

0-£© q Ul,...,U 9 e_R CFl=&q 

and the last sum over a is either equal to n x ei? n *( x )' or vanishes, depending 
on whether or not the u sequence is a permutation of the w sequence. Now the 
claim follows using (8). I 
We now apply Lemma 13 to get the following bound 

\l p \ < e -^> £ {x , y} ei l*-y| x K[ m] x Y[ (m(x) +m*(x))!3 

xe-R 




J CR d Vc[s(Z,h)} (^* ' VO V>0l) • ■ ■ ^(z, )V>* (Wi)...f (w ? ) 
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where we used the easily verifyable fact \m\ = \m*\ = i.e., in any J p one 
always has an equal number of ip and ip* factors remaining. Let us define the 
following initial multiplicities, i.e., before applying the derivatives, for the i/j 
and tp* fields respectively: 

mo(x) = \{ie = x and % = 0}| + 21 {xeT} + 2\{j, l<j< k\ yj = x}| 

and 

mS(x) = |{i € /|xi = x and Hi = *}l + 21 {xeT } + 2|{j, 1 < j < k\ Yj = x}| . 

Once again we have \mo\ = \thq\. When expanding the sums (32), each site will 
receive a number <5(x) of derivatives, as well as a number <5*(x) of 

derivatives. Besides, one trivially has 

£(x) + <f(x) = d(x) 

where d(x) denotes the degree of the vertex x £ -R in the tree T. 

Now it is easy to see that the maximal number of field factors occurs when 
all the derivatives pull new |^| 4 vertices from the exponential. Namely, for any 
x G R, 

m(x) < m (x) + 25* (x) + <5(x) , 
m*(x) < m^(x) + (5*(x) + 2(5(x) , 
m(x) + m*(x) < mo(x) + m^x) + 3d(x) . 

We will also need the identities 

^(m (x)+mS(x)) = |/|+4|T|+4fe (43) 



and 



which imply 



and 



d(x) = 2(|fl| - 1) (44) 



|m | + |mS| + 3|d| = |/| + 4|T| + 4k + 6|i?| - 6 (45) 



H < y + 2|T| + 2k + 3\R\ -3 . 



We are now ready to bound the number of derivation procedures p. 

Lemma 14 The number of derivation procedures is bounded by 

^2 1 < 2 l R l- 1 e l / l +4 l T l +4fc+6 l i? l- 3 J] d(x)! . 
P xe-R 
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Proof: One pays a factor 2' R I _1 for the sums (32). For x £ R\T the derivatives 
can only act on the fields already present, and we therefore have a number of 
choices limited by 

-n -n m (x)! rog(x)! 

J K'»oW>*W} Jl {"»SW>**W}( mo ( x )_ < y( x ))! ( m *( x ) - <5*(x))! ' 

For x G T, the number of terms produced by applying the derivatives is bounded 
by 

(46) 



d(x)! x exp <jm (x) + m^x) + -d(x) + 3 
Indeed, this is the number of terms produced when computing 



d 



5(x) 



9 



(x) 



5*(x) 



7/;(x) mo(x) V'*(x) mS(x) e -^(*)V(x) 2 



Let us first perform the 



x ) derivatives and then the 

a 



derivatives. When 



evaluating the very last derivative g^j^ one has to choose between deriving a 
new vertex from the exponential which gives a factor 2 for the choice of ip in the 
■0(x) 2 V'*(x) 2 vertex, or deriving a field factor which was already there which 
at most gives mo(x) + 25* (x) + <5(x) — 1 possibilities. Indeed, either the factor 
was present initially which corresponds to mo(x) choices, or it was in a vertex 
first derived from the exponential by a g^*^ derivative, in which case one has 
to pay a factor <5*(x) to identify that derivative and a factor 2 for the choice 
of field ip within the vertex — ^^(x.) 2 tp* (x) 2 . The last possibility is when the 
derived factor was in a vertex first produced by one of the previous <5(x) — 1 
derivatives of type ■ Such a derivative already consumes one of the two 
■ip's in the vertex, so we only have to pay a factor of <5(x) — 1. In sum, the last 
' derivative at most gives (mo(x) + 2<T(x) + (5(x) + l) possibilities. Likewise 

derivatives has at most (mo(x) + 25* (x) + <5(x)) options, 



the before last 



etc. Therefore, the number of possibilities for the 



di/>(x) 



derivatives is bounded 



by (m ?^^l 2 ( ?*^+i v 1)1 • B y a similar reasoning, that of the 



derivatives, 



(mo(x)+2<5*(x)+l)! • "J ^ lraB " nm B> ^ ^ av*(x) 

which are performed first, is bounded by ^^^j^yr-^ ■ As a result, the total 
number of possibilities is at most 

(m (x) + 25* (x) + <5(x) + 1) ! (m* (x) + 5* (x) + 1)! 



< 



(m (x) + 25*(x) + l)! 
m (x) + 2«5*(x) + d[X ' ] 1 ^ 



i(x) 



K(x) + i; 

<5*fx) 



mS(x) + 



3 

+ 2 



2 2 

where we used the arithmetic versus geometric mean inequality 
s! 



5*(x) 



(47) 



(s-q)\ 



= s(s-l) 
s + (s- 



< 



< 



■■(a -9 + 1) 
!) + ■■■ + (« 



+ D 



1 
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Finally using the inequality x n < n\e x for each of the two factors on the right 
hand side of (47), as well as the trivial inequality <5(x)!<5*(x)! < d(x)!, we find 
that the number of terms produced by the derivatives at the site x is bounded 
by 

d(x)! x exp|m (x) + m* (x) + ^<T(x) + ^<5(x) + 3 j . 

Now redo the same reasoning, but this time first applying the g^f^ and then 

the g^*^) derivatives. One will get the same bound but with <5(x) and <5*(x) 
exchanged. Taking the geometric mean of the two bounds gives the desired 
estimate (46). Using the same inequalities, the bound for x G R\T is easily 
seen to be no greater than the one for the x G T case. The lemma now follows 
from (43) and (44). ■ 
Putting the previous considerations together we now have a bound on X 
from (42): 

\X\ < 2l«l- 1 e l / l+ 4 l T l+ 4fc + 6 l«l- 3 ^ |/|+2|T|+2fc+3|i?h3 

xe -MoE {x , y}es |x-y| x -Q d ^y x -Q ( mo ( x ) +m *( x ) +3 d(x))!^ . 
xeR xeR 

Note that we used K i > 1 which is clear from (8) since X^ z ez d e _/1 ' z ' > 1. We 
now use the estimate 

(ra (x) +m^(x) +3d(x))! < (m (x) + m£(x))! x d(x)! 3 x 4™o(x)+m 5 (x)+3d(x) 
and identity (45) in order to write 

H (mo(x) +m*(x) + 3d(x))!l < 2 l / l+ 4 l T l+ 4fc + 6 l^- 6 

Yl d(x)!§ x J| (m (x) +m^(x))!^ . 



xeR 



We also need the following important remark: if \R\ > 2, then the integral X of 
(42) is zero unless \R\ < |T| + |7| + k. Indeed, if \R\ > 2 then the connecting 
tree T is nonempty and therefore each site x 6 R receives at least one field 
derivative. If such a site x contains no source, i.e., x / Xj, for any i £ I, and 
if the site is not in T, and if x / yj for any j, 1 < j < k, then the derivative 
has nothing to act on and the integral X is zero. Note that if \R\ = 1, then 
we already have a characteristic function in (41) enforcing |T| > 1 or |/| > 1. 
As a result, we always have \R\ < |T| + |7| + fc, which can be exploited by 
introducing the corresponding characteristic function. This condition, together 
with |T| > M — k and the assumption < A < 4 implies 

max{|_R|-|/|,M} 

(48) 



|T|+fc 



The previous considerations allow us to write 



d_ 

du 



M 



C(R,I,u\) 



< 1 



Vi G I 
xv G R 



Tci? 



\R\ > 2 or 
|T| > 1 or 
III > 1 
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max{|.R|-|7|,M} 



£ 1{|T| > M - A:} — 2 m U < |T| + | J| + A;} 

fc=o ' ^ ' 

Yl 2 ' Rhl x (2ev / ^i) |/|+4|T|+4fc+6|Rh3 x e" ME {-.y}6'^ |x " y| 

H d(x)!§ x (mo(x) + m$(x))!* (49) 



yi.-,y fc 6T 

'X tree 



xei? 



|T|! 



< 2l T l. 



where we used K\ > 1 as well as the inequality (M-k)\(\T\-M+k)l 

We now need a lemma which bounds local factorials of the degrees in the 
tree by a portion of the tree decay. This is a volume effect due to the finite 
dimensionality of the host lattice *Z d . 



Lemma 15 For any a > we have 



n <« 

xe-R 



x e 



|x-y| 



< K. 



\R\ 



where K 2 = maxj^.i, ^2,2} with 

K 2 ,i = 

and 



K 2 ,2 = exp J sup ^ 

rr6[l,+oo[ 



xlogx 



„ d ,d 
27T20(2 

r(i + i) 



«r (2 + 11 



2 2+ 



Q: 



d V7T 



Proof: We write the quantity to be estimated as the product over x G R of 

d(x)! e ~f ^yKx,y}«l x -yl (50) 

and we will bound the last expression using the fact that the d(x) sites y which 
are neighbors of x in the tree T are distinct, and the more they are the further 
away from x they have to be. Indeed, for any r > 0, the number B r of lattice 
points at distance at most r from x satisfies 



B r < 2° 



jz G N d , |z| < r}| 



< 2 d Vol J z+[0,l[ d 

\zSN d ,|z|<r y 

< 2 d Vol ({z G M^, |z| ^r + v^}) 
Vol(S d_1 ) y dp p**" 1 



< 



< 



d 

7T2 



r(f + i) 



(r + Vd) d . 
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Now if d(x) > 2B r , at least half of the d(x) neighbors of x are at a distance 
greater than r from x. This would imply 



e -fE y | { x, y }ei l*-yl < e -T d ( x ) 



Let us first suppose that 



Letting 



d(x) > 



27T2G(2 

r (i + !) 



r(f + 1) 



^ d 

2tt2 



d(xU -Vd, 



we see that r max > and r = r max satisfies the hypothesis d(x) > 2i? r . Conse- 
quently, it follows that 



e 2 Ey|{x, y }es l x y l < exp 



a:d(x) 



1 



^^d(x) -Vd 

27T2 / 



Using the trivial inequality d(x)! < d(x) d ( x ) and the definition of i^2,2 which 
clearly is finite when a > 0, we have 

d(x)! e ^f S y | { x, y}es l x -yl < K2 2 . 
In the second case where 

27T2d2 

d(x) < — - j 

r(| + i) 

then (50) is trivially bounded by i^2,i- H 
Since the d(x)! to be bounded in (49) appear at the power | and we want to 

use only a fraction, say half, of the tree decay; we will use the previous lemma 

with a = f so that fa = f. 

We now turn our attention to the (mo(x) + mj(x))! 2 in (49). Let 

c/(x) = |{i G 7|xi = x}| , 



6(x) 
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{xeT} , 
|{j,l<j<%,=x}| . 



Then 



[] (m (x) + m5(x))!5 = fj ( C/ (x) + 6(x) + 4u y (x))!* 

xe-R 

< [J { C/ (x)!6(x)!i; y (x)! 4 x 6 c/(x)+6(x)+4 % (x) j3 _ 



xeK 



xe-R 



That is 



! [ (mo(x) + m5(x))!3 < 6 i|/|+2|fl|+2k x2 4§l«l x "Q C/(x) ,| x TJ Uy(x) ,2 (51) 



xei? 



xeR 



xe_R 
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where we used |T| < \R\. 

We now simplify the bound (49) by using Lemma 15 with a = as well 
as (51). We also bound the sum over the no longer needed T by 2^1, and we 
bound k by M in the various exponents where it appears. After some cleaning 
up, we therefore get 



d 



M 



<( 4 



max{|_R|-|/|,M} 



yi £ i 

Xj G R 



M 



fc=0 

where 



xe-R 



£ n vyW 2 

yi,...,y fc e-R xei? 



V e -^E {x , y}es |x-y| 



% tree 



K 3 = 2 w 3 2 e 10 V6KfK* , 
K 4 = 2e v / 6^" , 
K 5 = 2 6 3 2 e 4 i\ 1 2 . 

Note that Ki in Lemma 15 is such that K<i > i^2,i > 1 and therefore K%, K4, K§ 
all are > 1. 

The sums over the positions yi, . . . ,y& of the vertices created by the ad- 
ditional Taylor expansion, as well as over the positions X2,...,x n in the l 1 - 
clustering bound will be done thanks to Lemma 11. 

Remark 6 We will use this lemma with both (5 = 2 and (5 = \ in order to sum 
over the yj and Xj respectively. However the (5 versus 1 dichotomy prevents 
us from compensating the factorials at the power 2 by the ones at the power \ . 
This is the main bottleneck we found on the way to Conjecture 1. 

Thanks to Lemma 11, the sum over the y's is bounded by 2\ R \ +k - 1 kl 2 . We 
then use the coarse bounds 

M 

k\ M\ 2l iJ l+ fe - 1 < M! 2 2l R l +M (M + 1) < Ml 2 2^+ 2M . 

k=0 

As a result we have the desired bound which is summarized in the following 
proposition. 

Proposition 3 Suppose the function J : Z d — > C satisfies Kq = 1, where Kq 
is the quantity defined in Lemma 3. Suppose M is a nonnegative integer, R is 
a polymer in A, u belongs to [0, 1] and < A < 4. Then, we have the estimate 



d 

du 



M 



((R,I,u\) 



Vi€l 
Xj G R 



< 



max{\R\-\I\,M} 



(2K 3 )^k{ II (4K 5 ) m M\ 2 



II C '( X ) ! " 



2 X 



xe_R 



e -imx,y}ed*-yl 



1 tree 



(52) 



where right derivatives are meant when u = 0, and left derivatives when u = 1. 
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4.2.2 The convergence criterion 



We now need to address the convergence criterion (38) of Proposition 2. Take 
any fixed site z G A, we then need to bound the quantity 

Q= Y HzeR}\C(RAm m ■ 
ReP(A) 

The simple case / = 0, u = 1, M = 0, of the basic raw estimate (52) gives 

Q< Y ^ zG ^}(^3) |R| Y e-^ E ^y}«l x -yl . 
ReP(A) |~>« 

% tree 

We now proceed as in §3.3 and condition the sum on m = \R\ and also introduce 
the identity (21), in order to eliminate R. The sum over the locations of the 
labeled sites gives a factor K™~ using (8) and letting K 6 = K x (d, The 
sum over trees t on the set [m], this time without the dt(*)!'s is simply bounded 
by m m ~ 2 < (m — 2)!e m by the coarser version of Cayley's Theorem [54, Prop. 
5.3.2]. Note that one has to treat separately m > 2 and m = 1. The result is 
easily seen to again be a geometric series bound, namely 

Q < Y (^K 3 K 6 ) m . 

m>l 

Therefore if A < K < 4, which we assume from now on, then we have 

1 



IIC(-,M)||< 2 



and the criterion is satisfied. 



4.2.3 Justification of the term by term differentiation 

We now address the issue of term by term differentiation leading to the series 
expression (40). Recall that as a corrolary of the mean value theorem and 
Lebesgue dominated convergence, the equation 

(1) N { ^ ^■■■^ tn ( x «))I«a = Yh E 4> T (Ri,...,R,) 

^ ' p>l Hi,...,Hp6P(A) 



6P(A) 



I q disjoint 



n 



E 



Ni+-+N P =N 



7Vi!...JV p ! 



n 



du 



((R q ,I q ,u\) (53) 



will be established for u G [0, 1] provided we can find majorants G(R, I, A, k) 
which are uniform in u such that 



d_ 

du 



< G(R,I,X,k) 
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for any u G [0,1], and such that for any integer M, 1 < M < N, one has 



E I E i#'(*....^)n{^ 

p>l y fl 1 ,...,fl p SP(A) 



E 



M! 



[jG(i? 9 ,/g,A,M g ) < +CX) . 



Mi!...M ! 

Ml+---+M P =M 1 p q=l 

We will take the majorants G(R, I, A, k) provided by the right-hand side of (52). 
We therefore have to show the finiteness of 



a = E^i E i* t w,...,k,p j / ' d,SJoml 

p>l & Si flp£P(A) 

/!,... ,/ p C[n] 



Ui, = [n] 



E 



Ml 



MiH \-M p =M 



M\\... ';/,! II 

F Q=l 



max{|fl,|-|/,|,Af,} 



(2^3)1^1^' (4if 5 ) M «M,! 2 



xi 



Vie/, 

X,; G iL 



n c/,(x)!* 



2 X 



xeK„ 



E 



|x-y| 



Since the only issue is finiteness, we will use very coarse bounds for Q. We 
write 

- a n max{|fl 9 |-|/ 9 |,Af,} / X -\\R q \-\h\ 

A) ~ U 

as well as 

v 

H H Cj ,(x)!2 <n!^ . 

9=1 XE-Rq 

We also bound the characteristic functions by 

J, disjoint 1 TT 1 / V< e M < 1 / 7 « r dis j°^ 1 1 {X1 € lRq } . 
Dl q = [n] J 11 \ Xi € R q / - \ U/, = [n] / 1 <>-i 5/ 



We finally use 



and 



and therefore get 



Mil . . . M p \ < Ml 2 M+P 

MiH \-M p =M 



(54) 



E 1 

A,...,ipC[n] 

Q < n! 



/.disjoint . _ 

u/, = [n] Sn - e 



2 ^ 4 M n M! 2 (8 ^ 5) M 



E 4 E e ^ =1 i?J|^ T (i? 1 , . . . , R p )\. P (Ri) ■ ■ ■ p(Rp) 



P>1 P " fii,...,R p 6P(A) 
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with 



% tree 



g-^F E {x ,y}e5 l x ^yl 



Now using the same argument as in §4.2.2 we see that 



m>l 



As a result, we will have \ \p\ \ < 1 provided 

1 



A < 



< 



1 



;i + 2e)eK 3 K 6 3eK 3 K { 



which we now assume. Finally, Lemma 6 shows that Q is finite, and therefore 
(53) as well as (40) are justified. 

4.2.4 The clustering estimate 

We now, under the hypotheses n > 2(iV + 1), Kq — 1 and < A < (i+2e)eK$K§ ' 
come to the clustering estimate proper, i.e., the bound on 



t= Y. H{ X 1 = 0} |(V ttl (xi),...,^"(Xn)) 

xi,...,x n eA 



T I 

a,aI 



Using (40), we can write 



(1-u) 



N-l 



xi,...,x„eA 



d " v nv-l)' ^ £ l^ 1 ^) 

p>i Ki,...,fl P eP(A) 



1 



iq disjoint 



E 



^ o=l 



A N v 

du 



C(Rq,I q ,u\) 



Ulq = [n] , 

q 1 1 ' Ni+-+N P =N 

Inserting the bound (52), performing the u integral, and tidying the resulting 
inequality, one obtains 

<C < (XK 5 fK2 Y ^=°}E^ E 



xi,...,x n eA 



pi 

P >1 1 i? 1 ,...,flpgP(A) 

'i 'pCH 



n 

9=1 



7Vi+-+JV p =JV 
>n max{|R 9 |-|/ 9 |-Af 9 ,0} ( w ■ _ r 



n «,(x)!3 



X < 



x£i? 



E 



-^E{x, y } e <5: 9 l*-y| 
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We now proceed as in §3.3 and push l{xi = 0} through the sums over p, the 
Rq 7 s and the I q 's, before bounding it by the coarser condition 1{0 € U^ =1 R q }. 
We likewise push the sums over the Xj's inside the appropriate bracket factor. 
The sums over the source localizations Xj are then performed using Lemma 11 
and yield the same bound as (25). Hence, 

£<(\K,) N (2K A TY;h E 1{0 e U^i?,} 



p. 

p>l 1 B lr ..,R p 6P(A) 
/l,...,JpC[n] 



l^i,...,^)!! 



I q disjoint 



^ )|/i|!...|/ P |!x Y, W- N J 



n 

9=1 



^\ ma,x{\R q \~\I q \-N q ,0} 



(4FT 3 ) x { 



N!+-+N P =N 
= -^E{x, y } e5 : g l x ~yl 



^ Tq tree 



We now introduce a constant 7, < 7 < 1, to be fine-tuned shortly; and we 
suppose that A satisfies the extra hypothesis < j < 7. We can now use the 
estimate 



mnx{\R q \-\I q \-N q ,0} 



< 7 max{|i? 9 |-|/ q |-7V 9 ,0} < 7 |fl,|-|/,|-iV, _ 



We also use the previously derived inequalities (27) and (54) which yield 



£ < 



N 



{^TN\n\ 

7 



E i E MO e u^ =1 i?J|^ T ( J R 1 , . . . ,/2p)| n ^) 



p. 

p>l ^ fl 1 ,...,fl p eP(A) 

Jl,...,JpC[n] 



9=1 



with 



ro(i?) = 4.(4 7 K 3 )I R I ^ e 



2 ^{x,y}e1 



|x-y| 



T tree 



Now, again as in §4.2.2, we have 

INI < 4 E (SnKsKef 



m>l 



Thus we will have II roll < 1 as soon as 



7< 



40eK 3 K 6 

We therefore choose 7 = 40ei ^^ 6 • Since we have by hypothesis that 



< A < 4 7 



< 



< 



10eK 3 K 6 (1 + 2e)eK 3 K 6 3eK 3 K 6 



< 4 
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we checked the validity of every condition we needed to check and we can use 
Lemma 6 and conclude 



{2\K 5 \ N ,4K 4 
7 



In other words we have the desired clustering bound 

£ |(^ 1 (0),V' fe (x2),...,^(x„))I| < 

x 2 ,...,x„eA 

A^ n\ N\ (80eK 3 K 5 K 6 ) N (160eK 3 K 4 K e ) n . 

We will now get rid of the restriction to Kq = 1 by a simple scaling trans- 
formation on the field variables V- Indeed if one does the change of variable 
ip = r]ijj f , for some rj > 0, in the original model, one easily sees that the n- 
point truncated correlations of the ip fields becomes rf 1 times the analogous 
correlation for the if/ fields. The latter are sampled according to the measure 
corresponding to the input parameters r] 2 J, 7] a \ instead of the original function 
J and coupling A respectively. Now one can go back through the definitions of 
the various constants fiQ, a, K±, K 2 , if 3, K 4 , K§ and Kq, and easily check that 
they are invariant if one multiplies the J function by a nonzero scalar rj 2 . On 
the other hand Kq = j ^f^y_j ^ gets multiplied by rf 2 . Therefore one can 

make the new Kq = 1 by choosing rj = j^(j(o)~j^) • Thus we have proved 
that if 

4(j(o) - J*? 

- 10ei^6</(0) 2 ' 
then, uniformly in A C Z d , and n even satisfying n > 2(N + 1), we have 

|(^(0), V' fc (x 2 ),...,^"(x n ))T| < Cl (N,J) xc 2 (J) n x\ N xn\ 

x2,...,x n eA 



where 



and 



Cl (N, J)=N\ 



( 80eK 3 K 5 K 6 J(O) 2 \ N 

{ 4( J (°) - J *) 2 ) 



C2 (J) = 160eK 3 K 4 K { 
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J(0) 



^(J(O)-J^) 



This completes the proof of Theorem 4, i.e., the clustering estimate for n-point 
functions with n > 4. H 
The proof of Theorem 5 is exactly the same as the one given above for 
Theorem 4 with the choices n = 2 and N = 1. The only difference is that the 
right-hand side of the starting point equation (40) now expresses the Taylor 
remainder 

(^(xi),^ 2 (x 2 ))I iA - (^( Xl ),^ 2 (x 2 ))I j0 
instead of the full 2-point function. 
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